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The main characteristics of the ground and, in particular, the first excited monopole 
state in the 4 He, 16 O and 40 Ca nuclei are studied within the generator coordinate 
method using Skyrme-type effective forces and three construction potentials, namely the 
harmonic-oscillator, the square-well and Woods-Saxon potentials. Calculations of den- 
sity distributions, radii, nucleon momentum distributions, natural orbitals, occupation 
, numbers and depletions of the Fermi sea, as well as of pair density and momentum dis- 

tributions are carried out. A comparison of these quantities for both ground and first 
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excited monopole states with the available empirical data and with the results of other 
theoretical methods are given and discussed in detail. 

PACS number(s): 21.60.-n, 21.10.-k, 21. 10. Re, 27.10.-h, 27.20. +n, 27.40. +z 



I. INTRODUCTION 

T— I ! 

VO 

The study of nuclcon-nuclcon correlation effects is important part of the contemporary nuclear physics 
[0J|]. The basic idea of the independent particle models (IPM) consists in the assumption that nucleons 
move independently in a mean field created by the same nucleons. Consequently, each coherent motion of 
the nucleons creates changes in the mean field. The Hartree-Fock approximation accounts only partially 
for the dynamic nucleon correlations. There exist many experimental data (see e.g. the review in |Q) 
showing that the IPM are unable to describe basic nuclear characteristics. For instance, the nucleon 
momentum and density distributions in nuclei cannot be reproduced simultaneously This is also 

the case with the occupation numbers, the hole-state spectral functions, with characteristics of nuclear 
reactions and others [0J|,|^-0|. This imposes the development of correlation methods going beyond the 
limits of the mean-field approximation (MFA) which account for nucleon-nucleon correlations (see e.g. 
[|| ) . This can be reached by extending the class of the trial functions which is used at the diagonalization 
of the nuclear Hamiltonian. 

The generator coordinate method (GCM) |^|| is one of the methods beyond the MFA which have been 
applied successfully to studies of the collective nuclear motions. For instance, the GCM has been applied 
to investigate the breathing-mode giant monopole resonance within the framework of the relativistic 
mean- field theory in jlCfl . In it the constrained incompressibility and the excitation energy of isoscalar 
giant monopole states were obtained for finite nuclei with various sets of Lagrangian parameters. An 
extension of the method of Ref. jl(| by using a more general ansatz for the generating functions of 
the GCM and by including the isovector giant monopole states has been done in [[ll]]. The use of the 
Skyrme effective forces has simplified the study of the monopole, dipole and quadrupole isoscalar and 
isovector vibrations in light double magic nuclei Jl2],[l3| . An approach to the GCM using square- well and 
harmonic-oscillator construction potentials has been applied to calculate the energies of the ground and 
first monopole excited state, as well as the density and the nucleon momentum distributions in the ground 
state of 4 He, 16 O and 40 Ca nuclei |l"4|]T5| ]. In [M the occupation numbers, the depletion of the Fermi sea 
and the natural orbitals (NO) which diagonalize the ground state one-body density matrix (ODM) have 
been calculated. The NO related to the ground state and the single-particle potentials corresponding to 
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them have been studied in detail in fl7fl . In [p|Jl8| the two-nucleon center-of-mass and relative motion 
momentum distributions of n — p pairs in the 4 He, 16 and 40 Ca nuclei have been calculated. The 
existence of high-momentum components in the one-nucleon and two-nucleon momentum distributions in 
the case of the square-well construction potential with infinite walls within the GCM has been obtained. 
The studies of the energies, the nucleon momentum and density distributions in 4 He and 16 nuclei have 
been extended by means of a two generator coordinate scheme in |i~9| ] . 

In the last years giant resonances with various multipolarities different from the well-known dipole 
resonance have been discovered. Among these collective excitations the monopole excitations take a 
particular place. The isoscalar giant monopole resonances or the so-called breathing vibrational states 
(with I n = + at energies of approximately 13 to 20 MeV) have been well established experimentally 
(e.g. [pp|-p7|) and have been considered to be compressional nuclear vibrations. Their study concerns 



the important problem of the compressibility of finite nuclei and nuclear matter (e.g. 1 28 - 33 1 ) . The 
description of such states is related mainly to general characteristics of the nuclei and weakly to the 
peculiarities of nuclear structure. Recently, in a series of articles Bishop et al.(see [Q and references 
therein) applied the translationally invariant cluster (TIC) method to light nuclei, in particular to the 
4 He nucleus. The basic properties, such as the energies and the density distributions of the ground and 
first excited breathing mode state in 4 He have been considered. 

The 4 He nucleus has a well-established spectrum of excited states. Calculations for the monopole 
oscillations of helium which practically involve the whole nuclear volume have been performed on the 
base of the nonlinear time-dependent Hartree-Fock (TDHF) approximation p5|. Its small amplitude 
limit, namely the random phase approximation (RPA), has been used extensively to describe nuclear 
collective motion |36|. The TDHF method |57| and the relativistic RPA J38| have been recently used to 
extend the study of isoscalar monopole modes in finite nuclei up to 208 Pb. It has been found that some 
effective Lagrangians can describe ground states and giant resonances as well, and in particular, they can 
predict correctly breathing mode energies in medium and heavy nuclei. In general, it has been pointed 
out that when going from heavy to lighter systems the trend is that the collectivity becomes weaker. This 
fact poses the question about the role of different kinds of nucleon-nucleon correlations corresponding to 
single-particle and collective motion modes and, therefore, the investigations on the correlation effects 
became an important task in the nuclear theory. 

The aim of the present work is to study the main characteristics of the ground and, in particular, 
of the first excited monopole state in the 4 He, 16 and 40 Ca nuclei within the GCM using Skyrme- 
type effective forces. It concerns the energies, the density distributions and radii, the nucleon momentum 
distributions, the natural orbitals and occupation numbers, as well as the pair center-of-mass and relative 
density and momentum distributions. As known, the natural orbital representation (NOR) ]39| ] gives a 
model-independent effective single-particle picture for the correlated states. We emphasize that in our 
work this is done also for the first excited monopole state in the nuclei considered on the basis of the 
ODM calculated for this state. It is known that the results of the GCM calculations depend on the 
type of the construction potential used to define the generating function in the method. In the present 
work we use three construction potentials, namely the harmonic oscillator ((HO) where the oscillator 
parameter is a generator coordinate), the square- well with infinite walls ((SW) where the radius of the 
well is a generator coordinate) and Woods-Saxon ((WS1) where the diffuseness of the Woods-Saxon well 
is a generator coordinate and (WS2) where the radial parameter of Woods-Saxon well is a generator 
coordinate). The values of the radial parameter in the WS1 case and the diffuseness parameter in the 
WS2 case are taken from pc[ . In our opinion, the calculations performed in this work can give an 
essential information about beyond MFA effects (accounted for in the considered approach to the GCM) 
on the mentioned characteristics. The results are compared with the available empirical data and with 
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the calculations of other theoretical methods. 

The basic relations of the GCM are given in Section 2. The results of the calculations and the discussion 
are presented in Section 3. The conclusions from the work are given in Section 4. 



II. BASIC RELATIONS IN GCM 

In the case of one real generator coordinate x the trial many-body wave function in the GCM is written 
as a linear combination || : 

/>oo 

*( ri )= / f(x)mr i },x)dx, i=l..A, (1) 

JQ 

where $({i"i}, a;) is the generating function, f(x) is the generator or weight function and A is the mass 
number of the nucleus. 

The application of the Ritz variational principle SE — leads to the Hill- Wheeler integral equation for 
the weight function: 

/>OC 

/ [H(x,x') - EI(x,x')]f(x')dx' = 0, (2) 
Jo 

where 

H(x,x') =< <S>({r l }x)\H\<S>({v 1 }x')>, (3) 

and 

I{x, x') =< $({rj x)\<S>{{vi} x') > (4) 

are the energy and overlap kernels, respectively, and H is the Hamiltonian of the system. 

If the generating function <&({ri} x) for a N = Z nucleus is a Slater determinant built up from single- 
particle wave functions <p^(r,x) corresponding to a given construction potential then the energy kernel 
(3) has the form @: 

H(x 1 x')=<<P({r. l }x)\<i>({r l }x')> J H{x, x', r)dr . (5) 

In the case of the Skyrme-like forces (for nucleus with Z = N and without Coulomb and spin-orbital 
interaction) H(x,x',r) is given by: 

H(x,x',r) = |-T+ ^ P 2 + -^(3*i + bt 2 ){ P T+f) + ^(9*i - 5t 2 )(X7p) 2 + ^t 3 p 2+ ° , (6) 

where to, t±, ti, ta, a are the Skyrme-force parameters. The density p, the kinetic energy T and the current 
density j are defined by: 

A/4 

p(x,x',r)=4 (^WaM^MO, (7) 

A/4 

T(x, x', r) = 4 (N- l )^V<pl(r, x) ■ V^(r, x') , (8) 

X,ju.=l 

A/A 

i(x, r) = 2 £ (N- 1 )^^, z)V^(r, x') - [V^(r, x)]^(r, x')} , (9) 
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where 

N x „(x,x') = J ^(r,s)^(r,a/)dr. (10) 
The overlap kernel (4) is given by: 

I{x,x') = [det(Nx^)} 4 . (11) 

Solving the Hill- Wheeler equation(2) one can obtain the solutions /o, /i, •■■ for the weight functions which 
correspond to the eigenvalues of the energy E ,Ei.... 

The one-body density matrix pi(r,r') of the ground (i = 0) and the first excited monopole (i = 1) 
states in this GCM scheme has the form: 



Pi (r, r') = J J fi(x)fi(xf)I{x, x')p(x, x',r, v')dxdx' , i = 0,l, (12) 

where 

A/4 

p{x,x',v,v')=AY j (N-^^l^x)^',*'). (13) 

It follows from (12) that the nucleon density distribution pi(r) and the nucleon momentum distribution 
rii(k) of the ground and the first excited monopole states can be expressed as: 

Pi(r) = J J fi(x)fi(x')I(x, x')p(x,x',r)dxdx' , i = 0,l (14) 

and 



where 



rij(k) = / / fi(x)fi(x')I(x,x')p(x,x' ,k)dxdx' , i = 0, 1 , (15) 



p(x,x',r) = p(x,x',r, r' = r) , (16) 



A/4 

p(x, x', k) = 4 (N- 1 )^^ x)^(k, x') (17) 



and <p(k, x) is the Fourier transform of ip(r, x). 

The rms radii can be calculated using the expression: 



r(i) = Jr* Pi (r)dr 
y J r z pi{r)dr 

As shown e.g. in |jl|,^2| the attractive properties of the single-particle description can be preserved in 
the correlation methods using the natural orbital representation |3^| . In it the one-body density matrix 
has the simple form: 

p,(r,r') = ^n«^ ) *(r)^ ) (r'), * = 0, 1 , (19) 

a 

where the natural orbitals (NO) ijja (r) form a complete orthonormal set of single-particle wave functions 
which diagonalize the density matrix pi(r,r') for the ground and the first excited monopole states. The 
natural occupation numbers (ON) ria for the state a satisfy the conditions: 
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0<nW<l, £n« = A. (20) 

a 

It is seen from (|l9|) that the natural orbitals and the occupation numbers can be found solving the 
equation: 

Pi {v,v')^{v')dv' = n^^\v), » = 0,1. (21) 

For nuclei with total spin J = the one-body density matrix can be diagonalized in the {I j m} subspace 
(l,j,m being the quantum numbers corresponding to the angular momentum, total momentum and its 
projection). In the case of nuclei with spherical symmetry the natural orbitals have the form: 

4$ m W=^(r)Y lm (0,p) = !£&Y lm (e,<p) ) i = 0,l, (22) 
^nlm( r ) f rom (22) and Pi(r, r') from (12) in equation (21) and the integration over the angular variables 



where R {(r) is the radial part of the NO's and Yi m (9,ip) is the spherical function. The substitution of 

&W ( 

nlm \ 

give the following equation for the radial part w£?(r) of the NO's and the occupation numbers : 

^ W (r,r>«(r')d/ = n$u<*>(r) , i = 0, 1 , (23) 

■jo 
where 

K® {r,r>) = JJ f^Mx'^x, x') (^'W'jK'llr. z)ft n *(r', , i = 0,l. (24) 

The summation in (24) is performed over all single-particle wave functions with given I forming the Slater 
determinant of the generating function and TZ n i(r, x) are the radial parts of these functions. 

The nucleon density distribution pi(r) and the nucleon momentum distribution rij(fc) for nuclei with 
J = can be written in the following form in the NOR: 



^w = iE 2 ( 2/+1 )E n Si^wi 2 ' ^m, (25) 

/ n 



47T 



where 

= (|) ia H)' rT 2 ji{kr)R^{r)dr, i = 0,l .27) 
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is the radial part of the NO in the momentum space and ji{kr) are spherical Bessel functions of order i. 
The depletion of the Fermi sea can be defined by the expression: 

Vil) = J £ (2i + l)(l-nW), i = 0,l, (28) 

where {.F} refers to the Fermi sea. 

The rms radius of the NO ip$ (r) is given by the expression: 



< r« >= 



fr 2 Itii^WlV , . 



The two-body density matrix p^> (£1, £2; £i, £2) m the GCM has the form: 
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P ( %i,6;£U2) = / dxf*(x) / dx'f(x')pW(x,x';t; 1 ,fr,t;' 1 ,&) , (30) 



where the coordinate £ includes the spatial coordinate r, as well as the spin (s) and isospin (r) variables. 
If the generating wave function is a Slater determinant constructed from a complete set of one-particle 
wave functions tpi(x,t;) the matrix p' 2 ' (x, x' '; £1, £ 2 ; £i, £2) can be expressed as [ j39[ : 

where 



< 2 W;£,o=£ (w-^^owo^eo, (32) 



p 



(JV 1 ) ;fc is the inverse matrix of: 

TV^x') =J2 f dv V * k (x^) Vl {x',0 , (33) 

and 

/(a;, a;') = det(A^) . (34) 

The two-particle emission experiments require some knowledge of physical quantities associated with the 
TDM. For example, using the two-body density matrix p^ in coordinate space, Eq.(p0|), one can define 
the pair center-of-mass local density distribution: 

P W (R) = J p( 2 ) (R + s/2, R - s/2)ds (35) 

and the pair relative local density distribution: 

p( 2 ) ( S ) = J (R + s/2, R - s/2)dR . (36) 

In momentum space the associated pair center-of-mass and relative momentum distributions can be 
defined: 

n (2) (K) = J n (2) (K/2+k,K/2-k)dk, (37) 

n W(k)=[ n^(K/2+k,K/2-k)dK. (38) 

The physical meaning of p^(s) and rt( 2 )(k) is the probability to find two particles displaced of a 
certain relative distance s = ri— r 2 or moving with relative momentum k = (kj^— k2)/2, respectively, while 
p*- 2 -* (R) and (K) represent the probability to find a pair of particles with center-of-mass coordinate 
R = (ri + r 2 ) /2 or center-of-mass momentum K = ki+k2. 

III. RESULTS AND DISCUSSION 

A. The ground and first excited collective state energies 

The Hill- Wheeler equation (2) is solved using a discretization procedure similar to that of Refs. |I^Jl3 |. 
The values of the Skyrme-force parameters in (6) in the case of square- well construction potential are the 
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same as in [^5|. They are determined to give an optimal fit to the binding energies of 4 He, 16 O and 40 Ca 
obtained from the Hill- Wheeler equation (2). The parameter set values to — — 2765, ii = 383.94, t% = 
—38.04, t 3 — 15865, a = 1/6 lead to the energies of the ground E and the first monopole excited E\ states 
(without the Coulomb energy) shown in Table |[ In the case of the harmonic-oscillator and Woods-Saxon 
construction potentials SkM* parameter set values (to = —2645, t\ = 410, <2 = —135, t% = 15595, a = 1/6) 
p3[ giving realistic binding energies are used. It can be seen from Table | that the energies obtained 
with WS1 are close to that obtained with HO, while the energies corresponding to WS2 are closer to that 
calculated with SW construction potential. 

The values of the excitation energies AE — E\ — Eq of the first monopole state (J 71 = + ) calculated 
within the GCM are given and compared with other calculations and some experimental data in Table 
||. It is seen that the GCM results with HO construction potential and in the WS1 case are in good 
agreement with the result of Brink and Nash Q] for 16 O nucleus and they are closer to the experimental 
values than the GCM results with the SW construction potential. The energy of the first excited + in 
4 He obtained in the present work is compared in Table [II] with the recent refined calculations of Bishop 
et al within the TIC method as well as with the results of the coupled rearrangement of the channels 
method @. 



B. The one-body density matrix p(r, r') of the ground and first excited monopole states 



1. The nucleon density distribution p(r) 



The one-body distribution function corresponding to the density distribution pi(r) from (14) is given 
by the expression: 

fli (r) =4^r 2 ( o 4 (r), i = 0, 1 . (39) 

The function gi(r) of the ground (i=0) and the first excited monopole (i=l) states of 4 He calculated in 
the GCM are compared in Figure 1 with those obtained in the TIC method. It can be seen that there is 
a large difference in the behaviour of the function g(r) for both states. This is due to the strong increase 



of the size of the nucleus, as is shown in Table III: the rms radius increases from 1.89 fm for the ground 
state to 3.00 fm for the first excited monopole state in the case of the HO, from 1.77 fm to 2.86 fm in 
the case of the SW and from 1.85 fm to 3.29 fm when using WS2. The significantly larger values of 
rms radii for the first excited state in A He nucleus in respect to their ground state values support the 
breathing-mode interpretation. At the same time, the nature of the first excited monopole state in A He 
is not still well understood. An attempt to solve this problem is made in Ref. ^(| where it is concluded 
that the i = 1 state is dominated by the 3 H + p configuration. The rms radii of the ground and of the 



first excited monopole states in 16 O and 40 Ca are also presented in Table IB. We note that the rms radii 
of the ground state in these nuclei obtained with Woods-Saxon construction potential are very close to 
the experimental values. 

For the first excited monopole state the one-body distribution function g\{r) calculated within the 
GCM using the WS2 potential and that one calculated within the TIC method are in good agreement. In 
contrast with this state, for a ground state both methods yield functions go (r) which differ significantly. 
This is due to the spurious centre-of-mass motion correction for 4 He which is accounted for in the TIC 
method J34|. At the same time effect doesn't influence substantially the first excited monopole-state 
density. 

The nucleon density distributions of 16 and 40 Ca of the ground po and the first excited monopole 
Pi states obtained within the GCM are presented in Figures 2 and 3, respectively. For the first excited 
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monopolc state the density distribution p\ (r) decreases in the central region of both nuclei and increases 
for large r. This behaviour corresponds to the increase of the rms radius in the first excited state. It can 
be seen also that the difference in the rms radii Ar rms = rfnL — rf'Xs decreases with the increase of the 
mass number. 



2. The nucleon momentum distribution n(k) 

The nucleon momentum distribution n(k) is one of the nuclear quantities which are sensitive to various 
types of nucleon-nucleon correlations. The momentum distributions for the ground and first excited 
monopole states of 16 have been calculated in the GCM using Eqs.(15) and (17) and are shown in 
Figure 4. As can be seen there is no essential difference between the behaviour of the n(k) in the ground 
and in the first excited monopole states. 

It is seen from Figure 4 that the nucleon momentum distribution n(k) obtained in the GCM depends on 
the construction potentials considered. The use of SW potential leads to an existence of high-momentum 
tail of n(k) which is not the case when the HO and WS potentials are used. Here we would like to note 
that: 

(i) the single-particle wave functions corresponding to the SW potential contain themselves high- 
momentum components due to the particular form of the potential and they are reflected in the GCM 
result, 

(ii) the lack of substantial high-momentum tail in n(k) when using HO and WS construction potential 
leads to the conclusion that the nucleon-nucleon correlations included in the GCM approach are not 
of short-range type. An additional study of the correlations accounted for in the GCM concerns some 
two-body characteristics and we will discuss the results for them in subsection C. 



3. The natural orbitals and occupation numbers 

The natural orbitals and occupation numbers have been calculated using Eqs.(21)-(24). A discretization 
procedure with respect to both r and r' has been applied solving Eq.(23) and the resulting matrix 
eigenvalue problem has been solved numerically. The natural orbitals obtained in the coordinate space 
are plotted in Figure 5 for the 16 nucleus. The corresponding rms radii and the natural occupation 
numbers of 4 He, 16 and 40 Ca are shown in Tables 4 and 5, respectively. 

One can see from Table |fy] and from the shape of the NO's given in Figure 5 that the rms radii of 
the natural hole-states increase from the ground state to the first excited state. The occupation numbers 
of the natural hole-states decrease in the first excited monopole state, while those corresponding to the 
natural particle-states increase (see Table 0). As can be seen from the same Table [vj the depletion of 
the Fermi sea increases in the first excited state. This is strongly expressed for the 4 He nucleus. The 
calculated depletion decreases with the increase of the mass number. 



C. The pair center-of-mass and relative density and momentum distributions 



The pair center-of-mass and relative density and momentum distributions Eqs.(35)-(p8|) are calculated 
by using of HO and SW potentials and are given in Figure 6 for the ground state of 16 nucleus. They 
are compared with the results obtained within the low-order approximation of the Jastrow correlation 
method (JCM) where the short-range correlations are explicitely involved. As it is seen from Figure 
6(a') the relative distributions p^(s) obtained in GCM and JCM have different behaviour in the region 
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from to 0.4 fm. This fact shows that the GCM doesn't take into account short-range but other kind of 
correlations which are obviously related with the collective motion of the nucleons. The curves of (K) 
calculated for both potentials decrease rapidly down at K>lfm _1 while the result obtained within the 
JCM shows a high-momentum tail. 

IV. CONCLUSION 

In the present work the characteristics of the ground and the first excited monopole states of the three 
nuclei 4 He, 16 and 40 Ca are studied within the GCM using different construction potentials. Though 
the results are sensitive to the type of the construction potentials used, some general conclusions can be 
summarized as follows 

(i) There is an increase of the rms radius in the first excited monopole state as a consequence of the 
increase of the density distribution for large r. 

(ii) There is not an essential difference between the behaviour of the nucleon momentum distribution 
in the ground and in the first excited monopole state. 

(iii) The use of the natural orbital representation makes it possible to study the effective single-particle 
picture not only of the ground, but also of the first excited monopole state. It is established that the 
depletion of the Fermi sea increases considerably in the first excited monopole state. This is strongly 
expressed in the case of 4 He. The calculated depletion decreases with the increase of the mass number. 

(iv) The study of the natural orbitals in the ground and the first excited monopole state shows that 
there are not substantial changes in their forms for a given nZ-state when they are calculated for the 
ground and the first excited monopole state. The rms radii of the hole-state natural orbitals are larger 
in the case of the first excited monopole state in comparison with those for the ground state. 

(v) The results on the one- and two-body density and momentum distributions, occupation probabil- 
ities and natural orbitals obtained within the GCM using various construction potentials show that the 
nucleon-nucleon correlations accounted in the approach are different from the short-range ones but are 
rather related to the collective motion of the nucleons. It turns out that the latter are also important in 
calculations of one- and two-body overlap functions which are necessary in the calculations of the cross- 
sections of one- and two- nucleon removal reactions. The work on the applications of the GCM overlap 
functions for the description of cross-section of (p,d), (e,e'p) and (7,p) reactions is in progress. 
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Figure Captions 



Figure 1. One-body distribution function (Eq.(p9|)) g(r) of 4 He calculated within the GCM using 
harmonic-oscillator (a), square- well (b) and Woods-Saxon (c) potentials and within the TIC method 



1 34 1 . The normalization is: J gi(r)dr = 1 (z = 0, 1) 



Figure 2. Nucleon density distributions of the ground (po) and the first excited monopole (p\) states 
of 16 O calculated within the GCM using harmonic-oscillator (a), square- well (b) and Woods-Saxon 
(c) potentials. The normalization is: J pi(r)dr — 16 (i = 0, 1). 

Figure 3. The same as in Figure 2 but for 40 Ca. 

Figure 4. Nucleon momentum distributions of the ground (no) and the first excited monopole (ni) states 
of 16 calculated within the GCM using harmonic-oscillator (a), square-well (b) and Woods-Saxon 
(c) potentials. The normalization is: J rii(k)k 2 dk = 16 (i = 0, 1). 

Figure 5. Natural orbitals for the ls-hole (a, a', a"), lp-hole (b, b', b") and 2s-particle (c, c', c") states 
of the ground (i — 0) and of the first excited monopole (i = 1) states of 16 O in coordinate space 
calculated within the GCM using harmonic-oscillator, square-well and Woods-Saxon potentials, 
respectively. 

Figure 6. The pair center-of-mass p^(R) (a) and relative p^^s) (a') density distributions and the 
pair center-of-mass n^ 2 ^(K) (6) and relative n^(k) (b') momentum distributions of the ground 
state of 16 calculated within the GCM using harmonic-oscillator and square-well potentials. The 
results obtained within the Jastrow correlation method ]48[ are also given. The normalization is: 
/ (K)dK = 1, / n( 2 ) (k)dk = 1, / (R)dR = A * {A - l)/2, / p^ (s)ds = A * {A - l)/2 
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TABLE I. Energies (in MeV) of the ground and the first excited monopole states in GCM with square-well, 
harmonic-oscillator and Woods-Saxon construction potentials (without the Coulomb energy). 



Nuclei 


Energies 


GCM 


GCM 


GCM 


GCM 








cw 

o vv 


vv Ol 


vv 


4 He 


E 


-29.51 


-37.10 


-26.31 


-28.03 




£1 


-8.06 


-9.87 


-5.38 


-8.29 


ie 


E 


-139.96 


-144.80 


-137.99 


-139.58 




Ei 


-115.01 


-111.66 


-114.09 


-112.50 


40 Ca 


E 


-410.07 


-404.23 


-409.23 


-411.40 




Ei 


-387.58 


-370.72 


-386.13 


-384.92 



TABLE II. 


The excitation energies of + breathing states (in MeV). 






4 He 




40 Ca 


GCM(HO) 


21.45 


24.95 


22.49 


GCM(SW) 


27.23 


33.14 


33.51 


GCM(WSl) 


20.93 


23.90 


23.10 


GCM(WS2) 


19.74 


27.08 


26.48 


Ref. | 




23.20 


27.30 


Ref. H 




24 




Ref. |o) 






18.3 


Ref. H ?7=8 




30.5 


29.5 


Ref. §§ 77=12 




30.0 


26.5 


Ref. §§ 


24.96 






Ref. H 


22.86 






Ref. @ 


26.49 


29.49 




EXP 




22.9 
(T=2) @ 


20.0 
(T=0) II 



TABLE III. Rms radii (in fm) of the ground and the first excited monopole states (first and second row, 
respectively) in 4 He, 16 and 40 Ca obtained with HO, SW, WS1 and WS2 construction potentials. 



Nuclei 


GCM 


GCM 


GCM 


GCM 




HO 


SW 


WS1 


WS2 


4 He 


1.89 


1.77 


2.28 


1.85 




3.00 


2.86 


3.43 


3.29 


16q 


2.67 


2.63 


2.70 


2.71 




2.87 


2.90 


2.89 


2.90 


40 Ca 


3.37 


3.40 


3.38 


3.39 




3.45 


3.52 


3.48 


3.51 
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TABLE IV. Rms radii (in fm) corresponding to the natural orbitals of the ground and the first excited monopole 
state calculated in the GCM with HO, SW, WSf and WS2 construction potentials. The Fermi level is denoted 

by (*)■ 



4 He-ground state 4 He-first excited monopole state 



state 


HO 


SW 


WS1 


WS2 


state 


HO 


SW 


WS1 


WS2 


3s 


3.44 


1.57 


5.45 


4.69 


3s 


3.40 


4.92 


2.68 


2.91 


2s 


1.78 


2.77 


2.62 


2.84 


2s 


2.80 


2.37 


3.03 


4.21 


Is* 


1.88 


1.76 


2.26 


1.84 


Is* 


3.03 


2.91 


3.44 


3.11 






ls O-ground state 






16 O- first excited monopole state 




state 


HO 


SW 


WS1 


WS2 


state 


HO 


SW 


WS1 


WS2 


2p 


3.32 


3.65 


3.82 


4.89 


2p 


3.79 


3.57 


4.11 


3.98 


2s 


2.97 


3.59 


2.68 


4.44 


2s 


3.36 


3.55 


3.82 


3.29 


lp* 


2.81 


2.70 


2.87 


2.85 


lp* 


2.97 


2.92 


2.97 


2.96 


Is 


2.18 


2.35 


2.13 


2.21 


Is 


2.29 


2.53 


2.19 


2.30 




40 Ca- ground state 






40 Ca- first excited monopole state 




state 


HO 


SW 


WS1 


WS2 


state 


HO 


SW 


WS1 


WS2 


3s 


4.22 


4.48 


4.92 


6.09 


3s 


4.57 


4.41 


5.03 


3.68 


2d 


4.28 


4.55 


4.88 


4.63 


2d 


4.58 


4.51 


5.00 


3.66 


2p 


3.88 


4.49 


4.28 


3.08 


2p 


4.14 


4.46 


4.86 


3.73 


2s* 


3.62 


3.56 


3.98 


3.85 


2s* 


3.71 


3.65 


4.02 


3.82 


Id 


3.64 


3.59 


3.63 


3.63 


Id 


3.70 


3.67 


3.71 


3.73 


lp 


3.08 


3.32 


3.06 


3.14 


lp 


3.13 


3.39 


3.09 


3.27 


Is 


2.42 


2.25 


2.04 


2.20 


Is 


2.40 


2.29 


2.13 


2.46 
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TABLE V. Occupation numbers and total depletion T> (Eq.(28)) of the ground and the first excited monopole 
state calculated in the GCM with HO, SW, WSf and WS2 construction potentials. The Fermi level is denoted 

by (*)■ 







4 He-ground state 






4 He-first excited monopole state 




state 


HO 


SW 


WSf 


WS2 


state 


HO 


SW 


WS1 


WS2 


3s 


< fO -3 


0.002 


< 10 -3 


< 10~ 3 


3s 


0.003 


0.009 


< 10" 3 


< 10" 3 


2s 


0.005 


0.014 


0.024 


0.006 


2s 


0.f49 


0.f52 


0.045 


0.142 


Is* 


0.995 


0.984 


0.975 


0.993 


Is* 


0.843 


0.886 


0.955 


0.8576 


V 


0.5% 


1.6% 


2.5% 


0.7% 


V 


f5.7% 


11.4% 


4.5% 


14.2% 






le O-ground state 






16 0-first excited monopole state 




state 


HO 


SW 


WS1 


WS2 


state 


HO 


SW 


WS1 


WS2 


2p 


< fO" 3 


0.013 


< 10" 3 


0.003 


2p 


0.067 


0.095 


0.08 


0.07 


2s 


< f0~ 3 


0.008 


< 10" 3 


0.003 


2s 


0.040 


0.057 


0.02 


0.07 


fp* 


0.999 


0.986 


1.000 


0.999 


lp* 


0.930 


0.903 


0.933 


0.926 


fs 


f.000 


0.992 


1.000 


0.999 


Is 


0.956 


0.942 


0.990 


0.955 


V 


< 0.1% 


1.2% 


< 0.1% 


< 0.1% 


V 


6.4% 


8.7% 


5.3% 


6.7% 






40 Ca-ground state 






40 Ca-first excited monopole state 




state 


HO 


SW 


WSf 


WS2 


state 


HO 


SW 


WS1 


WS2 


On 

OS 


< 1U 


0.014 


< 10 -3 


< 1U 


OS 


0.04f 


0.076 


0.038 


U.Uz4 


2d 


< 10~ 3 


O.OfO 


< 10~ 3 


< 10~ 3 


2d 


0.029 


0.055 


0.036 


0.024 


2p 


< f0~ 3 


0.007 


< 10" 3 


< 10~ 3 


2p 


0.02f 


0.039 


0.011 


0.026 


2s* 


0.999 


0.986 


1.000 


0.993 


2s* 


0.952 


0.923 


0.960 


0.969 


fd 


0.999 


0.990 


1.000 


0.993 


Id 


0.964 


0.945 


0.963 


0.970 


lp 


0.999 


0.993 


1.000 


0.993 


lp 


0.972 


0.961 


0.988 


0.968 


fs 


f.000 


f.000 


1.000 


0.994 


Is 


0.993 


1.000 


0.998 


0.994 


V 


< 0.1% 


0.9% 


< 0.1% 


0.6% 


V 


3.2% 


4.7% 


2.6% 


2.8% 
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